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1. Introduction 

Cellular Automata(CA) are simple mathematical models which can generate 
complex dynamical phenomena. In principle, CAs are spatially-extended discrete dy- 
namical systems which consist of a lattice of sites. The state of each site is in one of 
the g states at time t and has local interactions. The system evolve synchronously in 
discrete time steps according to identical rules. Cellular Automata was introduced by 
von Neumann in order to address the problems of self- reproduction and evolution 1 . 
However, in spite of their mathematical interests, CA has not attracted much at- 
tention until the introduction of John Conway's game of lifes 2 around 1970. About 
twenty years ago, Stephen Wolfram introduced the Cellular automata to the physics 
community as models of complex dynamical systems 3 and also as a new approach to 
the parallel computing scheme(For a review see Ref.4). The interest in CA's potential 
application continues to grow. In fact Cellular Automata have been used to simulate, 
for example, solar flares 5 , fluid dynamics 6 , crystal growth 7 , traffic flow 8 and galaxy 
formation 9 . 

One of the important question in CA studies is to determine the behavior of the 
system at large times. According to Wolfram's classification 3 the asymptotic behav- 
ior of CA can be classified into four classes which are : (1) spatially homogeneous 
state, (2) fixed or periodic temporal structure and (3) chaotic pattern throughout 
space. The fourth class of CA behaves in a much more complicated manner and was 
conjectured by Wolfram as capable for performing universal computation. 

Due to the complexity of cellular automata, numerical computations of time evo- 
lution become the main scheme in CA studies. In Wolfram's work 3 some statistical 
approaches was introduced to characterize quantitatively the patterns generated by 
CA evolution. Statistical methods such as entropies and dimensions can only provide 
average properties of cellular automata. However, exact results are always needed in 



2 



more elaborate discussions . To make this point more explicit let us consider the 
effects of noise on CA. Due to the fact that all physical systems are coupled to noisy 
environment, it is interesting to know if the above classification of the attractors 
of CA remains intact under the influence of noise. But this question can only be 
answered by knowing the deterministic CA exactly, and as a result exact evaluation 
of the attractor without noise is called for. Unfortunately such exact calculations only 
exists for a few cases. One of the most studied case is the GKL automaton proposed 
in the 80' 10 and surprisingly the detail analytic proof was completed about a decade 
later 11 . 

The analytical proof of GKL automaton constructed in Ref.ll was under the 
assumption that the system is infinite in size. However all numerical investigation is 
studied on finite lattice with periodic boundary condition. It is therefore important 
to investigate GKL model with periodic lattice. In this work we revisit this problem 
by modifying the method of Ref.ll and prove that the GKL automaton defined on 
a finite periodic lattice of size 2N can only have two uniform states if the initial 
configurations containing an uniform block of size larger than N. The best time 
estimation for reaching these uniform states is obtained. Our results also include 
previous proof special case. 

The plan for this paper is as follows. In the next section a brief review of the 
GKL automaton is given and several definitions are introduced. The main results of 
this work are given in terms of lemmas and theorem which are presented in section 
three. A brief discussion of our results are provided at the last section. 

2. The GKL Automaton on Finite Periodic Lattice 

The cellular automaton is defined on a finite size lattice with 2N sites. Any general 
configuration of CA will be denoted by a. The configuration of the system at time t 
is denoted by a* = {a f (x), 1 < x < 2N} where a l (x) which denotes the state on site 
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x takes on two possible values {+, — }. The periodic boundary condition is realized by 
a t (2N + k)= <7*(fc). The GKL rules are given by the following updating procedures: 




majority{a t (x), a\x + 1), a l (x + 3)} if o t {x) = + 
majority{a t (x), a f (x — 1), a f (x — 3)} z/ cr*(x) = — 



(1) 



Any configuration a can always be divided into two blocks which are denoted by 
ob and 07. ob is a block of size Nb , which contains cr(l) to o-(Nb), and 07 has 
size ATj, which contains cr(iV B + 1) to a(2N). In this work, N B is required to satisfy 
N B > N + 1 such that iV B > A7. The time evolution of a B is denoted by a B . In this 
work a block containing the same sign is called homogeneous. Furthermore, the initial 
configuration a always begins with a <tb which is a homogeneous contains only + 
signs whereas 07 can have arbitrary configuration . We will show in the next section 
that starting with this initial configuration the system will end up in an attractor 
which is homogeneous for the whole configuration with cr* = {cr(x) = +, 1 < x < 2N} 
in the limit of large t. That is to say that the block aj will be washed out at the 
end. Even though during time evolution the process can be quite complicate. This 
particular homogeneous configuration will be denoted by cr + . A similar attractor also 
exists for the case with initial condition containing a a B which is an homogeneous 
block of — sign. This attractor is denoted by a_ = {(t(x) = — , 1 < x < 2N}. Before 
we proceed it is useful to introduce some definitions which are in close resemblance 
to the ones employed in Ref.ll but with modification. 

Definition 1 a is good to a if it satisfies three conditions as follows: 

(1) For 1 < x < a — 2, if o~(x) = —, then a(x + 1) = +. 

(2) For 1 < x < a - A, if a(x) = -, then a[x + 3) = +. 

(3) Ifa(a - 3) = a(a) = -, then a (a - 1) = -. 
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With this definition, it is clear that any initial configuration, a is good to N B + 2. 
Furthermore, to the left side of a, the configuration a satisfying the above conditions 
contains neither two consecutive — nor two — at a distance 3, except near the bound- 
ary site a. That is, {+ H } is allowed if the rightmost site of this block is a. For 

example, the configuration a with 

o~ = {- + - + + + -- + + •••} 

is good to site 8, but not good to any site larger than 8. It is convenient to denote 

{ } as a{x; 2} with a(x) = a(x + 1) = — . This notation can be generalized 

as a{x; L} which denotes a block containing only — of size L where L > 2 and 
starting from <j(x) to a(x + L — 1). Besides, we also denote { — h H — } as a{x} with 
cr(x) = a(x + 3) = — . 

Definition 2 a is good if a is good for all sites. 

Definition 3 A block in a is perfectly good if this block contains neither a{x; L} 
nor a{x}. 

3. The Proof 

The main result in this section is to establish the fact that the system has two 
attractors which are the homogeneous states a + and <r_ respectively. First we show 
that perfectly good configuration a stays as perfectly good as the system evolves. 
This is state as a lemma: 

Lemma 1 // a 1 is perfectly good, then a t+1 is also perfectly good. 

Proof: It can be verified by direct checking that if a t+1 is not perfectly good, then 
a 1 is also not perfectly good. For example, if a t+1 is not perfectly good then it might 
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contain a{m;L}(A similar discussion can also be done for a{k}.). For this config- 
uration if cr*(m + 1) = +, then <j*(m) should be — such that either a l {m — 1) or 
cr*(m — 3) is equal to — . As a result, whatever the sign of cr*(m — 2) be, cr* must 
contain either o{m — 1; 2} or (x{m — 3}. The same argument can also apply for the 
case of cr*(m + 1) = — which implies the existence of either o{m — 1; 2} or <r{m — 2} 
in cr*. Therefore, for o t {m + 1) = ±, a* is also not perfectly good. 

It is useful to keep in mind that if o t {m) = — and a t+l (m) is also — , then a 1 
must contains either a{m — 1; 2} or cr{m — 3} on the left side of m. It is also inter- 
esting to know that the number of time steps for a B to retain perfectly good can be 
obtained and the result is proved in the following lemma. 

Lemma 2 For any initial configuration with an homogeneous a B , the perfectly good 
configuration of a B is retained during time evolution within N B — 2 time steps. 

Proof: It is noted that, if a^ 1 contains either a{m; L} or a{m}, then either a{m'; L'} 
or a{m'} are contained in o l B . The proof of this fact can be established in the same 
way as the proof of lemma 1 and is omitted. 

However there is one particular case which requires further discussion. As a l B x 
contains cx{3; L} or a{3}, there is a possibility that cr* contains a{2N}. If that 
happens then the proof is invalidate. To show that this particular case does not 
occur one has to establish that starting from any homogeneous a B , c B does not 
contain any <r{3; L} or <x{3} as t < N B — 2. 

It is obvious that, the homogeneous configuration a B can only be eroded from the 
right side at most one site at each time step. As a result, either <r{3; L} or 5"{3} may 
appear at t — N B — 2. For the case of a{3; L}, according to this at-most-one-site 
propagating property, the state of any site less than 4 must be + at t = N B — 3. 
Hence, a NB ' 3 (3) = + with a NB ' 3 (4) = a NB - 3 (Q) = -. However this can not lead to 
cr{3; L} at next time step. As a result, a NB ~ 2 does not contain a{3; L}. On the other 
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hand, for the existence of <t{3}, by using the propagating property, one can conclude 
in a similar fashion that ct Nb ~ 3 (3) = + and g Nb ~ z (A) = ct Nb - 3 (Q) = - which imply 
ct Nb ~ 3 (5) must be — . Such configuration is contradicted to the definition of a{3}. 
Therefore o Nb ~ 2 also does not contain a{3}. This completes the proof of this lemma. 

Combining the results of lemmas 1 and 2, one can conclude that within N B — 2 
time steps, if ob at time t is perfectly good, then ob at time t + 1 is also perfectly 
good. Moreover, Because our initial configuration of <jb is perfectly good, then in 
Nb — 2 time steps, ob retains its perfectly good configuration. We also mention that 
there is another way to see why lemma 2 only holds for t < N B — 2. As the — 
sign propagates to the site 3 at t — N B — 2, a NB ' 1 (3) is determined by the site of 
(t Nb ~ 2 {2N) which is contained in a - / and implies o~ Nb ~ 1 (2) = a NB ~ 1 (3) = — could 
occur. If that is the case then as is not perfectly good. 

Lemma 3 In N B — 2 time steps, if o l is good to a, then a t+1 is good to a + 1 . 

Proof: This lemma can be established by noting that if a t+1 is not good to a + 1, 
then o l is not good to a. For example, one can check that if any a{m; L} appears in 
a t+1 , then a t+1 is not good for all sites with x > m + 2. By using the same approach 
as one employed in proving the previous lemmas it can be shown that a 1 must contain 
either a{m — 1; 2} or a{m — 3}. This result implies a* is at least not good to m + 1. 
This conclusion can also be drawn for cr*(m + 1) = — . The case of cr{k} can also 
be treated in the same way and the conclusion remains the same as in the case of 
a{m;L}. According to our definition of perfectly good, lemma 2 ensures that a B 
keeps its perfectly good configuration in N B — 2 time-steps therefore Lemma 3 is only 
true for t < N B - 2. 

Lemma 4 If cr* is perfectly good and crl eed {m, L s } (with L s >3) is a homogeneous 
block of + sign in o l , then the length (L s ) of this block increases at least two sites at 
the next time step. 
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Proof: It can be proved directly. When a 1 is perfectly good, then it can contain any 
block which has all + signs with length other than 2. It is a trivial case for the block 
containing only one single + thus we just need to consider any block with the length 
larger than two. If o l contains a block cr* eed {m, L s }, then a t+1 contains a homoge- 
neous block of + sign starting from cr m (m) to a t+1 (m + L s — 2). Now consider the 
case with cr*(m — 1) = — , due to the restriction of perfectly good of a*, it is easy to 
check that cr*(ra — 2) and <7*(m — 4) must be +. Hence a t+l (m-l) = a m (m-2) = +. 
Furthermore, from lemma 1, a t+1 is also perfectly good therefore if o t [m — 3) = — , 
it must change sign with a t+1 {m — 3) = + at the next time step. However, the reader 
can easily check that there is no constraint on a t+1 (m — 4) and a t+1 {m + L 8 — 1). As 
a result, o"* eed {^; L s } becomes ultedi 771 — 3, L s + 2} in a t+1 . 

Equipping with these lemmas, one can now proceed to prove that the system 
will evolve into the attractor if the initial configuration has an homogeneous a B with 
N B = N+l. 

Theorem For a with a size of 2N, if the initial configuration has homogeneous 
<j b with N B >N+l then a N+Nl ~ 2 = a+. 

Proof: For our initial configuration of a B , <?° is good to N B + 2. According to 
lemma 3, a Nl ~ 2 is good to 2N . It is easy to check that the rightmost four sites 
of a N '~ 2 of all good configurations have + sign. These sites are on the left side 
of a B (since the lattice is periodic). Hence, in the next time step, one obtains 
a N i- 1 (2N — 1) = o Nl ~ l {2N — 2) = +. Furthermore, according to lemma 2, a B still 
maintains its perfectly good configuration at time Nj — 1 and therefore 07 contains 
neither a{m; L} nor cr{m}. Thus, a N '~ l is perfectly good. Moreover, at this time 
step, a B also contains a block a + with at least of size N B — (Nf — 1) > 3. As a result, 
there is at least a block a 1 ^^ 1 {!,?>} in a Nl ~ 1 with a perfectly good configuration. 



8 



From lemma 4, <7^,^ 1 {1, 3} will expand with the expansion rate of 2 sites per each 
time step. Therefore, the number of time step for cr^ 7-1 to become a + is at most 
2JV ~ 3+1 . As a result, with our particular initial configuration of a finite size lattice a, 
(J n i +n-2 _ ^ Thus the time for reaching the homogeneous state is obtained and is 
given by T < N r + N - 2. 

4. Discussions and Conclusions 

The results obtained in this work is that for any finite one dimensional lattice 
with periodic condition, the GKL automaton has only two homogeneous attractors 
if the initial configuration contained a homogeneous block of size bigger than half 
of the lattice. These results are closely related to the previous results which were 
obtained for infinite lattice 11 . The previous result studied the eroding effects which 
wash out any finite island embedded in any infinite homogeneous configuration. It is 
shown in this work that the infinite homogeneous configuration is not necessary. In 
fact the infinite homogeneous configuration can be replaced by a finite homogeneous 
block where elimination of the island also occurs as long as the homogeneous block is 
larger than half of the lattice. Furthermore the upper bound of the time for reaching 
the attractor is also obtained in this work. Even though the results obtained here is 
an improvement of the earlier results, however the story is not complete yet. This 
is due the fact that, on finite lattice, numerical simulations 4 has shown that for any 
initial configuration with a density of either + or — sign greater than 0.5 the system 
always ends up on the a± attractors respectively. The results in this work certainly 
confirm to this fact. Therefore it will be interesting to rigorously establish the results 
of simulation by relaxing the condition of a homogeneous Ob- Work along this line is 
now being investigated. 
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